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Abstract 

Recent measurements of the inelastic and total proton-proton cross section at the LHC, and 
at cosmic ray energies by the Auger experiment, have quantitatively confirmed fits to lower 
energy data constrained by the assumption that the proton is asymptotically a black disk of 
gluons. We show that data on p{p)p, Tr^p, and K^p forward scattering support the related 
expectation that the asymptotic behavior of all cross sections is flavor independent. By 
using the most recent measurements from ATLAS, CMS, TOTEM and Auger, we predict 
crfP (yi = 8 TeV) = 100.6 ± 2.9 mb and crf^ = 14 TeV) = 110.8 ± 3.5 mb, as weU as 
refine the total cross section a^^^{y/s = 57 TcV) = 139.6 ± 5.4 mb. Our analysis also predicts 
the total 7r^7r+ cross sections as a function of v^. 
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1 Introduction 



Recent high energy measurements of the inelastic proton-proton cross section, made possible 
by the Large Hadron Collider (LHC) and a new generation of cosmic ray experiments, have 
convincingly confirmed [1] indications [2, 3, 4, 5] in lower energy data that the total cross section 
(Ttot behaves asymptotically as the squared log of the center-of-mass energy y^, reminiscent of 
the energy dependence of Proissart's unitarity bound [6]. This energy dependence is now solidly 
anchored to all pp and pp total and inelastic cross section measurements, from threshold data 
averaged by finite energy sum rules, to the result at 57 TeV center-of-mass energy of the Auger 
cosmic ray array [5]. 

The energy dependence is suggestive of that predicted by an asymptotic black disk. Although 
the data itself does not cover asymptotic energies, from an extrapolation of the fits constrained 
by analyticity, the features of a black disk emerge, with a purely imaginary amplitude and a ratio 
of Oinel/^tot consistent with 0.5 within errors [1]. Additional, and independent, confirmation has 
been provided by LHC measurements of the shrinkage of the elastic scattering cross section 
[7]. From a parton point of view, the picture that emerges asymptotically is that of a proton 
composed of an increasing number of soft gluon constituents, each carrying a decreasing fraction 
of the proton energy. The asymptotic cross section, clearly emerging from available data is given 

by 

where M, historically identified with the mass of the pion, is now associated with the particles 
populating the Pomeron trajectory, i.e. glueballs. 

If one ascribes the origin of the asymptotic In^s term in pp and pp scattering to gluons only, 
then it is universal and its energy dependence as well as its normalization is the same for tttt, 
7rp, Kp, and 7p interactions via vector meson dominance. In other words, the role of quarks, 
and therefore the quantum numbers of hadrons, becomes negligible. Although there is still no 
rigorous derivation, the straightforward interpretation of the present data is that, asymptotically, 
particles of all flavors evolve into a universal black disk of gluons. The COMPETE Collaboration 
already proposed that this asymptotic behavior crtot — B log^(s/so) applies to all hadron total 
cross sections, with a universal value of the coefficient B [3, 8]. 

In order to empirically test this universality, the p{p)p, vr^p, and K^p forward scattering 
amplitudes are analyzed, and the values of B, denoted respectively as Bpp, B^^p, and Bxp, were 
estimated independently [9]. The analysis was refined [10] for Bxp- The resulting values are 
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consistent with the universahty, Bpp 2± B^^p 2± Bxp, and thus, the universahty of B is suggested. 
Recently strong indications for a universal and Proissart-like hadron-hadron total cross section 
at high energy are also obtained in the lattice QCD simulations [11]. In this work we first update 
the analysis of the p{p)p, vr^p, and K^p data by including newly measured LHC results as well 
as very high energy measurements based on cosmic ray data. Wc also fit the p{p)n data at the 
same time. Subsequently, assuming the universality of B, we calculate the tt^tt'^ total cross 
section cTtot^^is) at all energies. Similar analyses are also done in Refs. [12, 13, 14] by using 
different methods. 

Although challenging, the data on tt^tt"*" collisions could be extended to higher energies 
exploiting high intensity proton beam accelerator beams planned worldwide, such as Project 
X [15] of FNAL and J-PARC in Japan [16]. At a later stage these may develop into muon 
colliders. As an example. Project X, a high intensity proton source proposed at Fermilab, would 
deliver proton beams at energies ranging from 2.5 to 120 GeV [15] and secondary pion beams 
with -E(7r) « 2 — 15 GeV. A muon collider with Project-X-intensity pion beams would represent 
a ir~^TT~ collider with ^/s = 1 TeV and a luminosity of 10^^ cm~^/sec [17], not quite sufficient, 
even for measuring the large cross sections discussed here. Some manipulation of the secondary 
beams would be required. On the other hand, direct measurements of a^^^ in wide range of pion 
beam energy would be made possible. In the absence of such measurements we will extend our 
calculations of cr'^^'^^ (s) into the intermediate energy region using Regge theory. This will allow 
us to compare our predictions with indirect information [18, 19, 20, 12, 21, 22, 23, 24] extracted 
from processes such as 7r~p — >■ Tr~Tr'^n, 7r~7r~A++, assuming one-pion-exchange dominance. 

2 Update of the fits to (Jtotai 

2.1 Analysis of Forward p{p)p, tt^p, K^p,p{p)n Amplitudes 

The energy (momentum) of the beam in the laboratory system is denoted by ^{k). It is related 
to the center of mass energy by 

s = 2Mu + + m^, u = ^/k"^ + w? , (2) 

where m = M, fj,, uik for pp, np, Kp scattering, and M, fj,, and uik are proton, pion and kaon 
masses; respectively, s 2± 2Mu in high-energies. For p{p)n, M is replaced by neutron mass M„ 
and m = M. 

The crossing-even forward scattering amplitude, Fat\^)^ given by the sum of Pomeron and 
Reggeon (including P' trajectory) exchange terms, while the crossing-odd ^ {u) is given by 
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a single contribution from Reggeon (corresponding to vector-meson trajectories) exchange contri- 
butions. Here the subscripts ab and ab represent ab = pp, ir'^p, K'^p, pn and ab = pp, Tr^p, K~p, pn, 
respectively. We consider the exchange degenerate /2(1270)-, a2(1320)-trajectories for the 
crossing-even Reggeon (tensor-meson) term and the p-, cj-trajectories for the vector-meson term. 
Their imaginary parts are given explicitly by 

/™ Fi^V) = 4 + cf log^ + cS') + ^ , (3) 

V m m J m \mJ 

/™f.t'M = ^(-) . (4) 

where Cq', f3^, and /3y are unknown parameters in the Pomeron- Reggeon exchange model. The 
and cf' are introduced consistently with the Froissart bound to describe the increase of (Ttot 
at high energy. The intercepts are fixed with ax'(O) = 0.542, ayiO) = 0.455, which is taken to 
be the same as the Particle Data Group [8] . The amplitudes Im F^^^ (v) are related to the total 
cross sections ^"^^"''(s) by the optical theorem: 

,a6 



where a^^\s) ^ ^ Im Fj^^H^) ■ (5) 



In our analysis, p"'''"''(s), the ratios of real to imaginary parts of forward amplitudes, are 
fitted simultaneously with the data on a^oC'^. Real parts of the crossing-even/odd amplitudes 
are directly obtained from crossing symmetry F^^\e^'^i') = ±^^^^(1^)* as 

We introduce F^'^\o) as a subtraction constant in the dispersion relation [25]. The p"-'''°-^{s) are 
given by 

^a6,a6(g^ = Re F~"-^^"-\u) / Im F~"'^^°-\v), F~"-^'"-\v) = fI^\v) ± f'^~\v) . (8) 
2.2 Constrained Analysis with Universal Rise of (Ttot and Duality 

The contributions of the tensor term in Eq. (3) and the vector term of Eq. (4) to the ^"^^"^(s) 
are negligible in the high-energy limit -y/i — >■ 00, where they are well approximated by the d^-^ q 
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terms: 

<'t(s)-<ot(s) - B,,log^^ + Zab, (9) 

^0 

Air Att f pab 2\ 

where B,, = —cf, Z,, = — I c^' - ^) , (10) 

sf = 2Miyf + M'^ + m^, vf = me ^ , (11) 

where s'q is a scale for the coUision energy squared. By neglecting the small tensor-term con- 
tribution, <t"o^"* develops a minimum Z^b- B^h controls the increase of (t"q^"^(s) at high energy. 
In practise, tensor and vector contributions are negligible for -y/s ~ 50GeV where u^^^ and cr^o^ 
are described by Eq. (9). 

Two independent analyses [9, 10] of forward p{p)p, tt^p, K^p scattering using finite-energy 
sum rules (FESR) as constraints, demonstrated that the universality relation Bpp = Bt^p = Bxp 
is valid to within one standard deviation. In the present analysis, we include the p{p)n data and 
assume this universality from the beginning. 

Bpp = B-^p = Bxp = Bpn = B (12) 

It leads to constraints among c^^, c^^ , o^^ , and c^" from Eq. (10). 

Other powerful constraints are obtained from FESR [5] for crossing-even amplitudes, 

- j^^ ^^^F^t\^)d^ = ^ /_ <^^a^\k)dk = j_ {at,{k) + <Vfc))/2 dk , (13) 

where A^i^2 = -\/-^i,2^ ~ w^^- The integration limit A''2 is taken in the asymptotically high-energy 
region, while Ni is in the resonance-energy region. The left hand side of Eq. (13) is calculated 
analytically from the asymptotic formula of Im F^^^ given by Eq. (3) , while the right hand side 
is estimated from low-energy experimental data. Equation (13) imposes duality on the analysis. 
It allows us to constrain the high-energy asymptotic behavior with the very precise low energy 
data, through averaging of the resonances. 

Following Ref. [9], wc take A^i = 0.818, 5, 5 GeV for ab = Tip, pp, Kp, while N2 is commonly 
taken as iVa = 20 GeV. The FESR (13) yields [9] the constraints, 

(ttp) 102.2/3yP + 627.3c7 + 2572c7 + 10891c7 = 66.96 ± 0.04, (14) 
(Kp) 9.353/3^^ + 39.23c^^ + UA.lc^^ + 398.5cf ^ = 38.62 ± 0.07, (15) 
{pp) 3.481/3fr + 10.89cg^ -\- 27.50cf^ + 71.00c|^ = 90.38 ± 0.20. (16) 

The integrals of the experimental cross sections in the right hand side are estimated very accu- 
rately from low-energy data with errors less than 1%, and these equations can be regarded as 
exact constraints among parameters. 
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2.3 Updated Analysis Including LHC and Very High Energy Cosmic-Ray 
Data 

In order to determine the value of B more precisely, we now include three recent measurements, 
ATLAS, CMS and Auger, covering the very high-energy region in our fit: 

• ATLAS reported [26] app inelastic cross section crf^^^ at 7 TeV of 69.4ib2.4(exp.)ib6.9(extr.) 
mb where exp./extr. refers to errors from experimental/extrapolation uncertainties. By 
using the ratio <Ttot/o"mei at 7 TeV of 1.38, obtained from the eikonal model [27], af^^ 
is predicted to be CT^otC^ TeV) = 96.0 it 3.3 it 9.5 mb. Recently, this measurement was 
confirmed by the CMS collaboration [28] reporting (Jinei = 68.0 it 2.0(syst.) it 2.4(lum.) it 
4(extr.) mb, (where lum. refers to the error associated with the luminosity) giving crfg^. = 
94.0 lb 2.8 lb 3.3 ±5.5 mb at the same energy. We include these data omitting extrapolation 
errors. 

• The Auger [29] collaboration measured af^^^ at 57 TeV to be 90 ib 7(stat.) -ii (syst.) ib 
1.5(Glauber), where the last contribution to the error comes from Glauber theory. Using 

+12 

(^tot/(^inei = 1-45 at 57 TeV from Ref. [27], crf^t at 57 TeV is predicted to be 131 ±10 -16 ±2 
mb. We also include this result with statistical error only. 

• The TOTEM [30] has measured a total proton-proton cross section at -s/s = 7 TeV, 
98.3 ± 0.2(stat.) ± 2.8(syst.) mb. 

Experimental data of a^^f^ at A; > 20 GeV and p"^'"*' at A; > 5 GeV for p{p)p, ir^p, K^p 
scattering are analyzed. We also include the data of crf^'^" and ff^ at A; > 10 GeV. These data 
are fit simultaneously imposing on the parameters C2'i q-, P'^y ■, Fap\Q) the constraints (12) and 
(14-16). The highest energy data for a^C^ data reach 26.4(25.3) GeV for 7r-p(7r+p), 24.1 GeV 
for K'^p, 1.8 TeV for pp (Tevatron), 57 TeV for pp (Cosmic-Ray), and 23.0(26.4) GeV for pn{pp). 

Table 1: Best-fit parameters of the fit to dtot and /9-ratios of ir^p, K^p, and p{p)p scatterings. 
Constraints of the universality of B, Eq. (12) and FESR (14-16) are used. The brackets represent 
the most dominant uncertainties: the statistical errors for (3v and F^~^\0) and the systematic 
errors, which comes from the TOTEM measurement [30], for the other parameters. 



ab 


B{mh) 


Y^(GeV) 


^afo(mb) 


/3^' 






pp 
irp 
Kp 
pn 


0.293(26) 
0.293(26) 
0.293(26) 
0.293(26) 


4.64(88) 
5.10(73) 
5.18(76) 
12.00(75) 


34.63(65) 
20.72(39) 

17.76(43) 
38.90(26) 


6.44(35) 
0.143(12) 

0.408(95) 
2.67(34) 


4.393(41) 
0.0505(12) 
0.687(10) 
3.87(12) 


8.1(6) 
0.06(61) 
2.4(1.0) 
-15.6(6.8) 
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The number of parameters fit is 6 x 4 — 6 = 18. The fit is very successful despite the omission 
of systematic errors of the very high-energy data. The total is 'X^ /Ndf = 498.69/(604 — 18), 
with x^/Nd values of 225.25/245, 153.95/162, 63.84/111 and 55.64/86 iov p{p)p, tt^p, iiT^p, and 
p{p)n data, respectively. The results of our best fit to a^^f^ are shown in Fig. 1. The best-fit 

values of the parameters are given in Table 1. In order to estimate the systematic error of the 
universal value of B, we shift the central value of a^^^ at 7 TeV by TOTEM [30] as 98.3 it 2.8 mb. 
The corresponding variation of the best-fit value of B is regarded as the systematic error of B. 

B = 0.293 ± 0.004,tat ± 0.026,j,,t mb , (17) 

which is consistent with our previous estimates, 

B = 0.2817(64), 0.2792(59) mb [4] and B = 0.280(15) mb [9]. 

The systematic uncertainty of Eq. (17) is larger than the statistical error of our previous estimate. 
We consider this B value is the most conservative estimate from the present experimental data. 

3 The TTTT Total Cross Section 
3.1 Theoretical Predictions of cr^ot^^ 

We infer the (7^1'^^ (^) based on the analyses of forward p{p)p, Tr^p, and K^p scattering ampli- 
tudes. Based on the result of the previous section, we can predict a'^^'^^ at high energy. By 
using the relation s 2± 2Mz/ for ab = ap = pp,Kp,Trp in high-energies, (T^oi'^^i^) ^1- (^) 
be rewritten in the form 

afr"(5) = Blog'^ + Za, + P^'{^) ±Pv'{^] > (18) 

where 



So \sij \SlJ 



\ Si / 

si is introduced as a typical scale for the strong interactions which is taken to be si = 1 GeV^. 
It is natural to assume that the universality of B and sq extend to tttt scattering. The 7r^7r+ 
total cross sections a^^'^^ are expected to take the form 

/ \ aT(0)-l / \ av{0)-l 

<I'^+(s) = 5 1og2- + Z,, + /3fM- ±/3V(-] , (20) 

So \si/ \si/ 

where B and sq are given by Eq. (18). 
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The values of Z^^^^ Zxp-, Zpp in Table 1 approximately satisfy the ratios predicted by the 
quark model. 

Z^p : Zkp : Z^p = 20.72 : 17.76 : 34.63 - 2:2:3. (21) 

By using the quark model mcson/baryon ratio, Z-^t^ is Zj^-^ = ^Zj^p = 13.8 mb, while the 
Ztttt is also given by Z^tt = -^Zj^p = 12.4 mb, where the meson/baryon ratio is taken to 
be Z^^p/Zpp = 0.60 instead of 2/3. This assumes that the Z^b terms represent the conventional 
Pomeron exchange with a unit intercept (and no logarithmic terms) and that its coupling satisfies 
the Regge factorization. So our prediction is 

Z^^ = (12.4 ± 1.4) mb , (22) 

where the uncertainty is estimated from the difference between the above estimates. Actually 
this is the main source of uncertainty for our prediction at very high energy. Presently we have 
no rigorous theoretical way to estimate the accurate value of Sq'^, hence, we assume for simplicity 

= 5.10 ±0.73 GeV , (23) 

where the uncertainty comes from a difference between Sq^ and for our best fit given in 
Table 1. 

The coefficients f3^y take multiplicative forms in terms of Reggeon-aa(66) couplings jRaa,Rbb 
with = jTaa iTbb and /Jy* = 7v'aa Ivbb- In the case ah = pp and Kp, both /2(1270) 
and a2(1320)-trajectories contribute via the tensor-meson term and both p and w-trajectories 
contribute through the vector-meson term, while in the case ah = 7rp and tttt only the former 
trajectories contribute through the tensor and vector terms. For ab = pn, 02 and p contributions 
changes their signs from aa = pp case. Using Eq. (19), the values of are obtained from 
Table 1: 

^t' = 7/27r7r7/2PP = 19.4(1.6) mb, Py^ = Jpnnlppp = 6.11(14) mb, 

^T^= ^B^h,a2^RKKlRpp = 7.9(1.8) mb, PyP= Y.R=p,u^lRKKlRpp = 13.2(2) mb, ^24) 
~^T= Ihpp' + la^pp' = 46.4(2.5) mb, = ippp^^lu.pp' = 33.2(3) mb. 

Ihw^-la^w^ = 19.2(2.5) mb, = 7pp/ - 7a,pp' = 29.3(9) mb. 

The 7-couplings violates largely the relation of SU(2) fiavor symmetry: 7/2pp = 7a2ppi Ipw — 
luipp- Since 7p^7r/2 = ^ipKK = IujKK, 7/27r7r/2 = lf2KK = la2KK from SU(3), pljfy = 
is expected. However, it is also violated in Eq. (24). On the other hand, for vrvr scattering, 
Pt^ = 7|27r7r Py^ = 7p7r7r- They can be evaluated from the results oi p{p)p,p{p)n and ix^p 
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of Eq. (24): 



Pt — 



(19.4 ±1.6)^ 



(46.4(2.5) + 19.2(2.5))/2 
(6.11 ±0.14)2 



mb = 11.5 ± 0.9 mb, 



V 



mb = 19 lb 5mb , 



(25) 



(33.2(3) - 29.3(9))/2 

which are compared with the other estimates by using the same method applied to different 
inputs; (/3f^,;9^^) = (13.39,16.38) mb [31] and (8.95(24),21.8(9.0)) mb [14]. 



Table 2: Numerical values for the predictions of crj^t'^ ^^'^ their difference for a range of energies. 

Uncertainties from the errors of /3^\ decrease with the increasing energies, and become negligible 
above y/s ~ 40GcV, while the uncertainties from the errors on B(= 0.293 it 0.004st;at ±0.026syst 
mb) and ^/so{= 5.10 ± 0.73 GeV) become sizable above this energy. 



V^(GeV) 



3 
5 

10 

20 

40 

50 

100 

200 

500 

1000 



^rot" (mb) 



difference(mb) 



22.6±1.4^ ± 1.5 
18.3±1.4z±0.9 



I3t 



15.9±1.4z ± 0.4^^ 
16.0±1.4z ± O.4J 
18.1±1.4z±0.6so±0.4b 
19.1±1.4^ ibO.T^u ±0.65 
23.1±1.4z ±0.9^0 ± 0.9b 
28.3±1.4z ± 
37.1±1.4z ± 
45.1±1.4z ± 



11.2±1.4z ± 1.5^^ 
11.8±1.4z ±0.9^^ 
12.8±1.4z±0.4^^ 
14.6±1.4z ± 0.4 J 
17.4±1.4z ± 0.6,0 ± 0.4b 
18.6±1.4^±0.7,,„ ±0.6b 
22.8±1.4^ ± 0.930 ± 0.9s 

1.1.0 ± 
1.4,0 ± 2.2b 
1.6,0 ± 2.9b 



11.4±3.0^^ 
6.5±1.8^^ 



3.1±0.8^^ 
1.4±0.4^^ 



0.7±0.2 
0.5 
0.2 
0.0 
0.0 
0.0 



In summary, cr^'^'^^(s) are predicted by Eq. (20) with the parameters B from Eq. (17), 
Z^j^ from Eq. (22), s^^ from Eq. (23), and f3^y from Eq. (25). The numerical values of our 
predictions for several y^- values are given in Table 2. 

3.2 Comparison with Indirect Experiments 

There are no direct measurements of al^^^ at present, however, indirect data at low- and intermediate- 
energy have been extracted in Robertson73 [18], Biswas67 [19], Cohen73 [20], Pelaez03,04 
[12, 21], Zak;harov84 [22], Hanlon76 [23], Abramowicz80 [24]. They are compared with our 
prediction in Fig. 2. 

The /Nd values of our prediction with no free parameters for the whole data set with 
-y/s > 5 GeV are not good. We consider this comes from the quality of the data. Data are 
mutually in consistent even in the same collaboration using different method [22] . 

However, Our prediction for cr^t'^^' by considering the uncertainty from /3y^ shown by 
thin dashed line, seems to be consistent with the low-energy regions of Abramowicz80 [24] and 
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Zak;harov84 [22]. As was pointed in Refs. [12, 21, 14], these data have the natural connection 
to the low-energy data points by Pelaez03,04 [12, 21] which was recently updated in Ref. [32]. 

Our prediction is consistent with the recent other estimates [31, 14, 21]. The authors in 
Ref. [12] analyze cf^q^^ , <^tot''pp(~ (^tot + '^tot)/^) vr^vr^ data simultaneously. All the data, 
including [23, 24], as well as the data with very low energies, y/s = 1.38, 1.42 GeV [12, 21], 
are included in their fit. However, at energies of -y/s = 1.38, 1.42 GeV, the Regge theory is 
not guaranteed to work a priori, although it seems from Fig. 2 that it still provides a fairly 
good description at those low energies. This suggests that all sub-leading Regge effects, when 
combined, result in a rather small contribution. 

4 Discussion and Conclusion 

Our predictions for cr^t'^^ are shown along with the results of our best fit for cr^t^ ^"^^ 

Fig. 3. The difference in normalization of these curves is determined by the Zab-, -Zpp > Z^^p > 

Z^TT, while their increase with energy is described by the universal value of B, Eq. (17). 

There are a few comments as our concluding remarks: 

• We have previously predicted the a^^^ for LHC and cosmic-ray energies in [4, 9]. Now 
our previous predictions can be tested by using the new experimental data [26, 28, 29, 
30] shown in Table 3. The result of the fit in the present work is also shown, together 
with the predictions at y/s = 8, 14 TeV. Our predictions are in good agreement with the 
experiments. 

Our result is also compared with the other predictions [3, 13, 12] at s/s =14 TeV in Table 4. 
All the models give consistent results within their uncertainties. The central value of B in 

the present analysis becomes somewhat larger than our previous estimate as can be seen in 
Eq. (17). This larger B value comes from the TOTEM measurement [30], of which value 
includes a large systematic uncertainty. As a result our present prediction at y/s =14 TeV 
becomes almost the same as that of COMPETE collaboration as can be seen by Table 4. 
Our previous prediction based on duality constraint will be tested more strictly in the 
future LHC experiment s/s = 8 TeV. 

• The COMPETE collaboration assumed the universality of sq in their fit [3, 8]. We have 
tested this sq universality [9]. By applying the further constraints Sq^ = s^^^ = Sq^ in 
our analysis of p{p)p,Tr^p, K^p data(not including p{p)n data), we obtain B = 0.299(8) 
mb and the universal ^/sq = 5.59(30) GeV, which are consistent with B = 0.308(10) mb 
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Table 3: Comparison of our previous predictions of ^^^^(mb) with the new experiments by 
ATLAS [26], CMS [28], TOTEM [30] and Auger [29]. The result of the fit in the present work 
is also given. For the derivation of the experimental values, see the text. The numbers with 
parentheses are the fit results, and the others are predictions. 



V^(TeV) 


<t (BH) [4] 


<ot (11) [9] 


this work 


<ot(mb) (cxp.) 


7 


95.4±1.1 


96.0±1.4 


(98.2±2.7) 


96.0±3.3ea:p.±9.5ea:tr. : ATLAS 
M.0±2.8syst.±S.3Lum±^-5Extr. ■ CMS 

98.3±0.2,tot ± 2.8syst ■■ TOTEM 


8 


97.6±1.1 


98.2±1.5 


100.6±2.9 




10 


101.4±1.2 


102.0±1.7 


104.5±3.1 




14 


107.3±1.2 


108.0±1.9 


110.8±3.5 




57 


134.8±1.5 


135.5±3.1 


(139.6±5.4) 


+12 

131±105tat. -l6Syst.±^Glauber ■ Auger 



Table 4: Comparison of our prediction of a^^^{mh) at ^/s =14 TeV with the other works: 
COMPETE [3], AGGPSS [13], and FY [12]. * This uncertainty is quoted form the largest and 
smallest values given in Table I of Ref. [13]. 



V^(TeV) 


this work COMPETE [3] 


AGGPSS [13] 


PY [12] 


14 


110.8±3.5 111.5±1.2 -ii 


+10.2* 
100.3-12.5 


104±4, 113±4 



and ^ = 5.38(50) GeV of COMPETE collaboration [3, 8]. However, the value of 
this additional constraint on sq, x^/^df = 438.53/(517 — 11), becomes worse by 7 units 
(for extra 2 constraints) compared with our best fit X^/^DF = 431.48/(517 — 13) with no 
constraint on sq. The data seem to favor the fit without so-universality although the 
improvement is not remarkable in this case. 

• '^0^ ~ '^o" further assumed in the COMPETE analysis. In mini-jet modcl[33], the C2 
coefficient is described by gluon-gluon scattering, and thus it is the same as pp and pn 
system, while the ci coefficient includes the effect of quark-gluon scattering, thus, it is 
generically different between pp and pn. Correspondingly, Bpp = Bpn consistent with the 
universality, while ^ Sq". Our best fit value \/io" ~ -'-^-^ GeV is slightly larger than 
"v/sq^ = 4.6 GeV. If we take Sq" = as one more constraint, the resulting for the 
[p^pn data is [Ndf = 69.3/(86 — 4) which is 14 units larger than the original value 
X^(pn data), x^/^df = 55.6/(86 — 5). The reduced x^ is less than unity also for the 
So-universal fit, but the experimental data prefer the non SQ-universal fit. So we did not 
adopt the so-universality in the present analysis. 
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otot Imb] 




Figure 1: Fit to the data of a^^^ (blue triangles) and a^^^ (black circles). The solid lines are 
our best fit and the dashed lines correspond to one standard deviation of B. The vertical line 
on the X-axis represents the lowest energy of the fit region y/s > 6.27 GeV corresponding to 
A; > 20 GeV. The LHC ATLAS [26] and CMS [28] data (with no extrapolation errors) at 7 
TeV (orange) and the Auger data [29] (with only its statistical error) at 57 TeV(orange). The 
TOTEM [30] at 7 TeV is shown by red. 
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Figure 2: vr^vr"^ total cross section (mb) versus -y/s. cr^t'^ (thick solid green) and cr^t'^ (thick 
solid black). Thin dashed lines represent the uncertainty from Py^ , which is the largest in the 
relevant energy region. 
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Figure 3: Total cross sections(mb) versus ^/s. Solid lines are pp,pp, iT~p, tt~^p, tt~ 7r~^ , it~^ 7r~^ from 
up-to-down. pp,pp and 7r~p, vr+p are our best fit, while 7r~7r+, tt'^'tt"'" are our predictions. The 
dot-dashed lines represent the upper(lower)-limit of our predictions of cr^^t^ {f^tot^ )• Dashed 
lines for p{p)p and n^p represent the uncertainties of our predictions which are obtained from 
the errors of B and Z^tt. 
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